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Abstract — The codeword stabilized ("CWS") quantum codes 
formalism presents a unifying approach to both additive 
and nonadditive quantum error-correcting codes (arXiv:quant- 
ph/0708.1021). This formalism reduces the problem of construct- 
ing such quantum codes to finding a binary classical code 
correcting an error pattern induced by a graph state. Finding 
such a classical code can be very difficult. Here, we consider 
an algorithm which maps the search for CWS codes to a 
problem of identifying maximum cliques in a graph. While 
solving this problem is in general very hard, we provide three 
structure theorems which reduce the search space, specifying 
certain admissible and optimal ((n,K,d)) additive codes. In 
particular, we find there does not exist any ((7,3,3)) CWS 
code though the linear programing bound does not rule it out. 
The complexity of the CWS-search algorithm is compared with 
the contrasting method introduced by Aggarwal & Calderbank 
(arXiv:cs/0610159). 



I. Introduction 

Quantum error correcting codes play a significant role 
in quantum computation and quantum information. While 
considerable understanding has now been obtained for a broad 
class of quantum codes, almost all of this has focused on 
stabilizer codes, the quantum analogues of classical additive 
codes. Recently, a number of nonadditive quantum codes have 
been discovered, with superior coding parameters ((n, K, d)), 
the number of physical qubits being n, the dimension of the 
encoded space K, and the code distance d [1], [2], [3]. These 
new codes have inspired a search for more high-performance 
non-additive quantum codes [4], a desire to understand how 
non-additive codes relate to additive codes, and how these may 
be understood through a cohesive set of basic principles. 

A systematic construction, providing a unifying approach to 
both additive and nonadditive quantum error-correcting codes, 
has been obtained [1]. This codeword stabilized quantum 
codes ("CWS" quantum codes) approach constructs the de- 
sired quantum code based on a binary classical code C, chosen 
to correct a certain error pattern induced by a self-dual additive 
quantum code which is without loss of generality, taken to be 
a graph state Q. The construction thus reduces the problem of 
finding a quantum code into a problem of finding a certain 
classical code. All previously known nonadditive codes [5], 
[6], [2], [7] with good parameters can be constructed within 
the CWS construction. 

The natural challenge in these approaches is efficient iden- 
tification of suitable classical codes, from which the desired 
additive and non-additive quantum codes can be constructed. 
It is apparent that due to the error pattern induced by the 
graph state Q, the binary classical code C does not coincide 
with the usual binary classical code where the minimum 



Hamming distance is a more important code parameter - 
although interestingly, they do coincide in the special case 
where Q is an unconnected graph, so the family of CWS 
quantum codes includes classical ("bit-flip") codes as depicted 
in Fig. 1. 




Fig. 1. The relationship of CWS codes with additive quantum codes and 
classical codes: ALL: all quantum codes; CWS: CWS codes; ADD: additive 
codes; CLA: classical codes. 

The CWS construction, observing that a classical code 
correcting certain bit-flip error patterns gives rise to a quantum 
code, allows a natural encoding of the problem of finding 
a quantum code Q = (G,C) into an equivalent problem, of 
finding the maximum clique of an induced graph, called the 
CWS clique graph. The existence of such a mapping is not 
surprising, since MAXCLIQUE is an NP-complete problem [8], 
[9], and thus can be used for a reduction from all unstructured 
search problems. In practice, many heuristic and randomized 
Clique solvers and SAT solvers have been developed, with 
reasonable run-times for small problem sizes. And since the 
search for CWS codes starts from a graph state Q, prior art 
in categorizing local Clifford (LC) orbits of those states [10], 
[11] helps simplify the problem. Nevertheless, without further 
simplification, a mapping of the CWS quantum codes search 
problem to MAXCLIQUE leaves the problem unsolved, due to 
the exponential computational cost of solving MAXCLIQUE. 
The real situation is even worse. For a general graph state, the 
search problem is NP-complete due to the reduction to MAX- 
CLIQUE. However, to search for all the quantum codes, we 
need to search for all graphs of n vertices, which contributes 
a factor of order 2" . 

Here, we present an algorithm for finding CWS codes, based 
on a mapping to MAXCLIQUE. We show that despite the expo- 
nential complexity of solving this CWS-MAXCLIQUE problem, 
the algorithm can be usefully employed to locate and identify 
a wide variety of codes, by taking careful steps to prune 
the search space. In particular, we show how the complexity 
cost can be reduced by using known graph isomorphisms and 
LC equivalences of graph states. We also present simplifying 
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criteria for the search, arising from the structural properties 
of CWS codes. We prove three theorems limiting whether 
((n, K, d)) additive codes with optimal K can be improved, 
or not, by the CWS construction. These theorems allow 
significant practical reduction of the search space involved 
in finding CWS codes using CWS-MAXCLIQUE. Furthermore, 
these theorems also indicate the existence of quantum codes 
outside of the CWS construction, as alluded to in Fig. 1. 

We also compare and contrast the CWS codes with another 
framework ("AC06") which was introduced independently [12] 
and is based on a correspondence between Boolean functions 
and projection operators. We interpret the AC06 framework 
to use a quantum state and a classical code, to generate the 
desired quantum code, but in a sense, it works in the reverse 
direction, starting from the classical code and obtaining the 
quantum state. We show how the AC06 Boolean function / is 
the analogue of our classical code C, up to a LC equivalence. 
This allows us to extend AC06 to degenerate codes, and to 
show that the AC06 framework can also be used to construct 
a search algorithm for new quantum codes, with comparable 
complexity to cws-maxclique. 

II. The cws-maxclique algorithm 

The CWS-MAXCLIQUE algorithm is a procedure to search 
for a quantum error correction code Q = (Q, C), given a graph 
state Q which maps quantum errors £ in the Pauli group into 
binary error patterns, and a classical code C, which corrects the 
error patterns. We present this algorithm below, beginning with 
a review of the basic definitions of CWS codes, proceeding 
to the details of the procedure, then rounding up with an 
evaluation of the computational complexity of the algorithm. 

A. Non-degenerate and degenerate CWS codes 

The basic concepts and definitions of CWS codes are 
described in a previous paper[l], and may be summarized as 
follows. The standard form CWS code is fully characterized 
by a graph Q and a classical binary code C, such that 
the corresponding CWS code may be denoted by the pair 
Q = (G,C). We define 

Clg{£) ={Cl g {E) I Ee£} (1) 

as the set of classical errors induced by quantum errors £ 
acting on the graph Q; these are the errors that the classical 
code C must detect. For each quantum error E, it is sufficient 
to express E in Pauli form as E = ±Z V X U for some bit 
strings u and v. The mapping to classical error strings is 

n 

Clg(E = ±Z v X u )=v®($(u) l r l , (2) 

where r; is the Zth row of the adjacency matrix for Q, and 
(u); is the I th bit of u. 

Using these definitions, the main theorem of the CWS code 
construction (Theorem 3 of [1]) may be given as: 

Theorem 1: A standard form CWS code, Q = (G,C) for 
graph state Q and classical code C, detects errors from £ if 



and only if C detects errors from Clg{£) and in addition, for 
each E E £, 

either Clg(E) ^ (3) 
or Vi Z C 'E = EZ C ' , (4) 

where Z Ci are codeword operators for C from {Z c } c& c- 

The case where Clg(E) ^ for all E £ £ is the non- 
degenerate case. For degenerate CWS codes, it will be useful 
to introduce a new set of classical bitstrings 

D g {£) = {c £ {0, 1}" | Clg(E) = and (5) 
c • u ^ for some E = ±Z V X U £ £} . (6) 

These bitstrings indicate codewords which are inadmissible, 
because they violate the condition given by equations (3) and 
(4) of Theorem 1 . Specifically, fix a codeword c, then for all 
£e£we must have Z C E = EZ C if Clg(E) = 0. Writing 
E = ±Z V X U , c is not an admissible codeword if Clg(E) = 
and c • u 7^ 0. In other words, if a CWS code is degenerate, 
some low weight errors act trivially on the code space (i.e. 
Clg(E) = 0), and these errors must act trivially on each 
basis state generated from the graph state Q (i.e. [Z c , E] = 0). 
Dg(£) describes basis states for which this is not the case. 

B. The CWS-MAXCLIQUE algorithm 

Given a graph Q, the problem of finding a CWS code Q = 
(Q,C), which corrects for quantum errors £, is reduced to a 
search for suitable classical codes. It is thus natural to ask 
how such classical codes can be found. One solution might be 
to use existing classical codes for this construction. However, 
that approach gives sub-optimal code parameters, due to the 
fact that C should be able to detect errors of the highest weight 
of the induced error patterns in Clg{£). This means that the 
classical code C must have distance significantly greater than 
that of the corresponding quantum code (G,C), as shown in 
the following example: 

Example 1: Let Q be an n qubit ring graph. If £ is the set 
of single qubit Pauli X, Y , and Z errors, then the induced 
classical errors Clg{£) are single, triple, and double bit flips 
respectively. Choosing the classical code C to be a binary 
((n, K, 7)) code results in a CWS code (Q, C) with parameters 
((n, K, 3)). However, C also detects many additional errors 
which are unnecessary for this construction, such as all the 
one to six bit flip errors; Clg(£) only includes a subset of 
those errors. 

This example motivates a search for specific classical codes 
which correct just the relevent errors for the CWS construc- 
tion. However, classical coding theory provides no efficient, 
systematic constructions for codes that correct the potentially 
exotic error patterns involved in the CWS construction. On the 
other hand, finding a code with the best K for given n and 
d is a problem which can be naturally encoded into an NP- 
complete problem such as MAXCLIQUE. This classic approach 
has been employed, for example, to show that the (10, if, 3) 
classical code with A' = 72 has optimal parameters [13]. 

CWS-MAXCLIQUE is a mapping onto MAXCLIQUE, of the 
problem of finding the CWS code (G,C) with the largest 
possible dimension K, for given parameters n, d, and graph 
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Q. The CWS-MAXCLIQUE algorithm gives steps to solve this 
problem, and is given in detail in the Algorithm 3 box. It 
proceeds in several simple steps. The first step, Setup (5, A) 
(Algorithm 1), finds the elements of Clg(£) and Dg(£). The 
second step, MakeCWSCliqueGraph(CL, D) (Algorithm 2), 
constructs a graph, denoted as the CWS "clique graph," whose 
vertices are classical codewords and whose edges indicate 
codewords that can be in the same classical code together. 
When searching for ordinary classical codes using an analo- 
gous procedure, the usual condition for joining two vertices 
by an edge is that the vertices are Hamming distance d apart. 
In our situation, vertices are joined by an edge if there is no 
error induced by the graph state that maps one codeword to the 
other. Finally, an external subroutine findMaxClique(V, E) 
is called; this routine is to employ known techniques to find 
the maximum clique in the CWS clique graph. The clique- 
finding subroutine is not specified here because there are many 
exact and heuristic techniques known in the community, for 
solving this classic NP-complete problem. Note that in the 
detailed description of the algorithms, two functions are used: 
String (i) : integer i — > binary string of i with length n, 
and its inverse, Integer(z) : binary string with length n i ^ 
integer of i. Also, an error configuration is a list of ordered 
pairs (LOC,TYPE) where LOC is the coordinate of the 
affected qubit and TYPE is one of X, Y, or Z. 

C. The complexity 

CWS-MAXCLIQUE is not an efficient algorithm; the run-time 
is at least of order ~ 2", because of the representation of the 
bit-string sets Clg(£) and Dg(£). These are needed to specify 
the CWS clique graph, which has 2™ nodes. In principle, 
instead of storing all this in memory, the vertices and edges of 
this graph could be computed on the fly, during execution of 
the flndMaxClique subroutine. However, these inefficiencies 
are not limiting factors, because of the even larger size of the 
search space involved in typical applications. 

Typically, the goal is not to search for an optimal CWS 
code, given Q and £, but rather, to determine if an ((n, K, d)) 
code exists when n and K are fixed. When K is fixed, finding 
a maximum clique is not necessary; rather, a clique of size 
K is desired. There are ( K ) such possible cliques. Checking 
whether a size K subgraph of a CWS clique graph is a clique 
just requires checking if that subgraph is fully connected. 
Given an adjacency matrix for the CWS clique graph (and 
constant time access to the matrix elements), checking a 
subgraph takes order K 2 steps. 

Searching over the space of all possible graphs Q involves 
searching a space of graphs with n vertices, with a total of 
possibilities. Therefore, the complexity of searching for 
an ((n,K,d)) CWS code is roughly 



(7) 



However, several practical improvements allow this search 
space to be pruned usefully. First, not all graphs Q need 
be considered; only those which are inequivalent under local 
Clifford (LC) operations need be checked. The LC orbits 



Algorithm 1 Setup(£,A): Compute Clg{£) and Dg(£), 
where £ is a set of Pauli errors and A is the adjacency matrix 
associated with graph Q. 

Require: A T = A, A^ = {0, 1} and A u = 
Ensure: CL[t] = <J(String(«) G Cl g {£)) and D[i] = 
5(StringW G D g {£)) 

i: for i G {0, 1}" do 

2: CL[Integer(i)] <- 

3: D[Integer(i)] «- 

4: end for 

5: for error configuration E G £ do 

6: ERR^- String(O) 

7: ERRX< — String(O) 

8: for (loc, type) in E do 

9: if type is X or Y then 

10: ERR <— ERR © (row LOC of A) 

11: ERRX <— ERR © String(2 L0C ) 

12: end if 

13: if type is Z or Y then 

14: ERR <— ERR © String(2 L0C ) 

15: end if 

16: end for 

17: CL[Integer(ERR)] «— 1 

18: if Integer(ERR) is then 

19: for i G {0, 1}™ do 

20: if errx ■ i ^0 then 

21: D[i] <- 1 

22: end if 

23: end for 

24: end if 

25: end for 

26: return (CL,D) 



Algorithm 2 MakeCWSCliqueGraph(CL, D): Construct a 
graph whose vertices V are classical codewords and whose 
edges E connect codewords that can belong to the same 
classical code, according to the error model indicated by 
Clg{£) and Dg{£). 



Require: CL and D are binary arrays of length 2" 
Ensure: 0" eV, n ^ v eV ^ D[v] = and CL[v] 
(v, w) G E CL[v © w] = 

1: V <— {0™} 

2: E <- 

3: for s G {0, 1}™ do 

4: if D[s] = and CL[s] = then 

5: V^VU{s} 

6: for v G V \ {s} do 

7: if CL[v © s] = then 



= 0, 



E 

end if 
end for 
end if 
end for 
return (V,E) 



EU{(v,s)} 
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Algorithm 3 CWS-MAXCLIQUE (£, A): Find a quantum 
code Q detecting errors in £, and providing the largest possible 
dimension K for the given input. The input A specifies the 
adjacency matrix of the graph Q. The output C is a classical 
code such that Q = {G,C) is a CWS code detecting errors in 

£. 

Require: A T = A, A tj = {0, 1} and A,, = Vi 
Ensure: K = \C\ is as large as possible for the given input, 
0" G C, and C satisfies the conditions in Theorem 3 of [1] 

l: (CL,D) <- Setup (£, A) 

2: (V, E) <- MakeCWSCliqueGraph(CL, D) 

3: C <- findMaxClique(F,£;) 

4: return C 



of graphs are well understood, and efficient algorithms exist 
to check for LC equivalence [10], [11], [14]. Therefore, the 
factor 2( 2 ) can be significantly reduced. A lower bound on 
the number of LC inequivalent graphs is given in [15], based 
on the number of non-isomorphic tree graphs, which roughly 
scales as 3". This reduction has played a key role in allowing 
us to employ the cws-maxclique algorithm on spaces 
with parameters up to n = 11 and K = 32. However, no 
suitable upper bound is presently known, which would give a 
quantitative estimate of the extent of the search space reduction 
due to LC equivalence. 

A second practical improvement comes from intrinsic prop- 
erties of CWS codes, which rule out existence of codes of 
certain ((n,K,d)) parameters, and relate the existence of 
certain parameter values with the existence of others. We will 
return to discuss these structure theorems in Section IV. 

III. Boolean functions and Classical Codes 

The CWS construction unifies all known additive and non- 
additive quantum error correction codes of good parameters, 
including both degenerate and non-degenerate codes. An al- 
ternative framework ("AC06") for non-degenerate codes, has 
been presented by Aggarwal & Calderbank [12], based on 
a correspondence between Boolean functions and projection 
operators. Because AC06 implies a search algorithm for quan- 
tum codes which is in a sense the reverse of that employed 
above, in CWS-MAXCLIQUE, it is interesting to consider the 
differences. 

In this section we study the relationship between AC06 
and the CWS construction, by linking the AC06 Boolean 
function, which we interpret to specify a certain classical 
code, to the classical code C used in the CWS construction. 
The components of the AC06 construction can be naturally 
associated with those of the CWS construction. In this way, 
we show that AC06 codes are spanned by a set of stabilizer 
states generated from a single state and a set of Pauli operators. 
Therefore, AC06 codes can be described completely, and in 
our opinion more transparently, as CWS codes. 

That this identification between AC06 and CWS is natural 
was mentioned previously [1], but the transform required has 
not been presented before. It is well known that any stabilizer 
state is equivalent under some LC transform to a graph state. 



Thus, supposing that a local Clifford operation maps the AC06 
stabilizer state to a graph state, it would be nice if this Clifford 
also described a transform from the Boolean function / to the 
binary classical code C of the CWS construction. Below, we 
show this mapping indeed exists, up to a technical subtlety 
with regard to the choice of the generating set for the stabilizer. 

The AC06 framework is not entirely complete since de- 
generate codes cannot be described as presented in [12]. 
Degenerate codes may, in some cases, outperform the best 
known nondegenerate codes. Such an example may be pro- 
vided by the [[25,1,9]] code obtained by concatenating the 
[[5, 1, 3]] code, since this is the best known [[25, 1]] code, it is 
degenerate, there is no known nondegenerate [[25, 1,9]], and 
it has the highest possible minimum distance [16]. We take 
the constraints given for degenerate codes in the CWS con- 
struction and map these backwards to given new constraints 
for degenerate codes in the AC06 framework. 

Given a complete AC06 framework which includes both 
non-degenerate and degenerate codes, we can then compare 
and contrast the computational cost of the CWS and AC06 
approaches for seeking optimal parameter quantum codes. 
When the search goal is to find an optimal ((n, K, d)) code 
for fixed n and K, the AC06 framework seems at first 
to involve a search over possibly 2 2 Boolean functions, 
while CWS-MAXCLIQUE involves a search over 2» possible 
graphs. This appears to give significant advantage to CWS- 
MAXCLIQUE. However, we find that with careful analysis of 
AC06, and extending it include degenerate codes, the two 
search algorithms have comparable complexity. 

A. AC06 quantum error-correcting codes are CWS codes 

A 11- variable Boolean function is a mapping / : {0, 1}" — > 
{0, 1} that maps a binary n-vector v = (vi, . . . , v n ) to a bit 
f(vi, . . . ,v n ). A Boolean function is nonzero if there exists 
some v such that /(v) = 1. We know that a Boolean function 
is naturally associated with a classical code 

C/ = {cG{0,l} n I /(c) = 1}. (8) 

A nonzero Boolean function / can be represented as 

/(v)= £ (9) 
cec, 

where v\ = Vi and = Vi = Vi © 1. The summation is taken 
to be modulo 2, i.e. XOR. The weight of a Boolean function 
/ ^ \C f \. 

The complementary set of a nonzero n-variable Boolean 
function /(v) is defined by 

Csetf = {a G {0, l}' 1 I X)/(c)/(c9a)=0}. (10) 

CGC7 

We know that the complementarly set is simply the set of 
vectors a such that Cf D (Cf © a) = 0, i.e. it is the set 
of (classical) detectable errors of Cf, since no codeword is 
mapped back into the code by a. 

Definition 1 (Definition 6 of [12]): Let P and Q be projec- 
tion operators on a Hilbert space H with K = image (P) and 
L = image(Q). Then 
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. P < Q iff K c L and K ^ L 

• P V Q is the projection of H onto the span K V L, the 
smallest subspace of H containing both K and L 

• P A Q is the projection of H onto K P\ L 

• Pis the projection of P onto K- 1 
. P©Q = (PAQ)V(PAQ). 

Definition 2 (Definition 7 of [12]): Given an arbitrary 
Boolean function f(vi,...,v n ), the projection function 
/(Pi, P2, . • • , P n ) is the expression in which Vi in the 
Boolean function is replaced by the projection operator Pi, 
multiplication (AND) in the Boolean logic is replaced by the 
meet operation P V Q in the projection logic, summation 
(OR) in the Boolean logic is replaced by the join operation 
P A Q in the projection logic, and the NOT operation in 
the Boolean logic is replaced by the not operation P in 
the projection logic. Note that summation modulo 2 (XOR) 
is replaced by the cooresponding operation P Q in the 
projection logic. 

Theorem 2 (Theorem 1 of [12]): If (Pj_, P 2 , . . . , P n ) are 
pairwise commutative projection operators of dimension 2 n_1 
such that (PxP 2 ...P„), (PP2...P0, (P 1 P 2 ...P n ) 
are all one-dimensional projection operators and H is of 
dimension 2", then Pf = /(Pi, P 2 , . . . , P n ) is an orthogonal 
projection on a subspace of dimension K = Tr(P/) = wt(/). 

Let (a|b) denote the concatenation of two ?i-bit binary 
vectors a and b. The symplectic inner product of 2n-bit binary 
vectors (a|b) and (a'|b') is 



(a|b)©(a'|b') = (a|b) 



/ 
/ 

= a • b' © a' • b. 



(a'|b') 



(ID 

(12) 



The symplectic weight of a vector (a|b) is the number of 
indices i at which either ai or bi is nonzero. P( a |t>) is defined 
by ei <E> e 2 <£> • • • ® e n where a equals / if (flj, 6j) = (0, 0), X 
if (a,, h) = (1,0), Z if (ai,bi) = (0,1), and Y if (a 4 ,6,) = 
(1, 1) and the associated projector is P( a |b) = h(I + P(a|b))- 
The next definition specifies the ingredients of an AC06 
quantum error-correcting code (AC06 QECC). Theorem 1 of 
[12] defines a quantum code, but our definition of an AC06 
QECC is based instead on Theorem 2 of [12], which provides 
sufficient conditions for the code to be an error-correcting 
code. 

Definition 3 (AC06 QECC): Let / be an n variable 
Boolean function and let x\, x 2 , ■ ■ ■ , X 2n be a list of the 
77-bit column vectors of an n x 2n matrix Af. An AC06 
QECC with data (/, {^i}f=i) is the image of the projector 
/(Pi, P2, . . . , P n ), where (i) the rows of Af are linearly 
independent with pairwise symplectic inner product zero and 
(ii) Pi = P( ai |bi) is associated to the 7th row of Af. 

Theorem 3 (Theorem 2 of [12]): Let be the set of all 
2n-bit vectors of symplectic weight less than d. An AC06 
QECC with data (/i{^t}|"i) is an ((n,K,d)) quantum code 
if / has weight K and {AfW T | w £ Dd} C Csetf. 

The main result of this subsection, stated and proven next, 
is that AC06 QECCs are CWS codes. 

Theorem 4: An AC06 quantum error-correcting code is a 
codeword stabilized quantum code. 



Proof: Consider an AC06 QECC with data (/, {x^i)- 
The matrix Af, whose 2n columns are {xi}?™i> nas linearly 
independent rows with pairwise symplectic inner products 
that are zero. Therefore, Af corresponds naturally to a group 
generated by n pairwise commuting operators {gi}™ =1 from 
the 77 qubit Pauli group. Let \S C ) be the state stabilized 
by S = ((— l) Ci 5i)" = i for some 77-bit vector c. A nonzero 
Boolean function / can be represented as 



cec f 



which corresponds, in this case, to the projector 

f(p u p 2 ,...,p n )= J2 p?p?---p, 



(13) 



(14) 



where P? = P t = \{I - g. t ) and P/ = P = i(J + 9i ). 
The term P^P^ 2 . . . P° n projects onto the state \S e ), where 
c = c~ic~2 . . . c n , therefore 



f(P 1 ,P 2 ,...,P n )= E \S C ){S C \. 



(15) 



cec. 



Hence, the AC06 QECC is spanned by a set of eigenstates of a 
stabilizer S, each of which has a vector of eigenvalues given by 
a codeword b in the inverted code Cf, where bi = indicates a 
+1 eigenvalue for gi and bi = 1 indicates a —1 eigenvalue for 
gi. To establish correspondence with a CWS code, we need to 
show that there is a mapping W from n-bit strings c to Pauli 
operators W(c) such that \S C ) = W(c)\Sqo...q). Indeed, there 
is a Clifford circuit U that encodes U\ 00 . . . ) = |Soo...o) and 

n 

acts like UZiU^ = gi for i = l,...,n. Therefore, UXiU^ 
anticommutes with gi and commutes with all gj, j ^ i. By 
this observation, the map 



W(c) :=\{[UX t U^ 



(16) 



has the desired properties, and we obtain the set of CWS 
word operators W(Cf) by applying W to each codeword 
in Cf. Therefore, the AC06 QECC with data {f,{ x i}i=i) 
is associated with a CWS code (not in standard form) with 
stabilizer state IS") corresponding to Af, classical code Cf, 
and word operators W(Cf). ■ 

The mapping can be inverted to obtain data for an AC06 
QECC from a CWS code as well. There is freedom in the 
choice of generating set for the stabilizer state in the CWS 
construction so it may be necessary to conjugate by a Pauli 
operator to fix the signs of the stabilizer generators to +1 
before mapping them to the column vectors {xi}^ v 

Example 2: This detailed example demonstrates the map- 
ping given in the proof of Theorem 4 from an AC06 QECC 
(/,{^}f=i) = {LAf) to a CWS code (S A ,C, W(C f )). The 
AC06 ((5, 6, 2)) code is given by the boolean function 

f(v) = V1V2V3 + V 3 V4V 5 + V 2 V 3 V 4 

+ ViV 2 V 5 + ViV4V 5 + V 2 V 3 V4V 5 
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and the matrix 



1 1 1 1 1 



1 

1 
1 

1000111011 
1 1 1 



110 111 
11000111 



First, consider the boolean function /. Indeed, f(v) is a 
function of n = 5 variables and has weight K = 6. This 
can be seen by writing / in the form 

f(v)= J2 /(cK 1 •••<"= 

c£{0,l}'* c6C/ 

= VxV 2 v s v A v 5 + Viv 2 v 3 v 4 v 5 + v~iv 2 v 3 v 4 v 5 

+ V1V2V3V4V5 + VlV 2 V 3 V 4 V 5 + V~lV 2 V 3 V 4 V 5 

where equals v,i if Cj = 1 and Vi if Cj = 0. The classical 
code Cf is the set of n-bit strings on which / evaluates to 1, 
i.e. 11100, 00111, OHIO, 11001, 10011, and 01111. Second, 
observe that the rows of Af are indeed linearly independent 
and pairwise orthogonal in the symplectic inner product. The 
rows of Af correspond to stabilizer generators E\ = IZYYZ, 
E 2 = ZYYZI, E 3 = YYZIZ, E 4 = YZIZY, and 
E$ = IZIXX, respectively. These are the generators of the 
stabilizer Sa for the state IS"). The AC06 construction uses 
the fact that the projectors P y = ^(1 + E y ), y = 1, . . . ,n, 
are pairwise commutative projection operators of dimension 
2"" 1 and P 1 P 2 . . . P n , P 1 P 2 ...P n ,..., PiP 2 ...P n are all 1- 
dimensional projection operators, so that Pf := /(Pi, . . . , P n ) 
is a projector onto a subspace of dimension wt(f) (Theorem 
1 of [12]), where the boolean operations are replaced by the 
operations defined in Definition 6 of [12]. Considering just the 
first term of Pf, we see that 

Pi f\P 2 A P3 A P4 A P5 

= PiP 2 P 3 {I - P 4 )(/ - P 5 ) 

= ^(1+ Ei)(I + E 2 ){I + E 3 )(I - E 4 )(I - Es) 

is a projector onto a stabilizer state where Wi is 

a Pauli operator that commutes with {Ei,E 2 ,E 3 } and 
anticommutes with {E 4 ,E§}, i.e. Wi = Z5. Notice that the 
partition of the generators into commuting and anticommuting 
sets is given by the first codeword 11100 of Cf. The terms are 
combined using the operation P(BQ = P + Q — 2PQ, which 
equals P + Q when the projectors are pairwise orthogonal, 
as they are when P and Q project onto stabilizer states. 
Therefore, P f = Y,f=i Wi\S){S\W} where the W, are 
chosen to commute or anticommute with the generators of 
the stabilizer of 15) according to the codewords of Cf. We 
conclude that the AC06 ((5,6,2)) code is a CWS code with 
stabilizer {IZYYZ, ZYYZI, YYZIZ, YZIZY, IZIXX) 
and word operators {Z§, Z 3 , Z 4 , Zi, Z 2 , X 3 X 4 X5 } 
that correspond to the classical code C — Cf = 
{00011, 11000, 10001, 00110, 01100, 10000} specifying 
the generator's signs for each basis state of the quantum 
code. We can arrange for the all-zeros codeword to be in 
C by multiplying each word operator by X 3 X 4 X$ (and, 
hence, adding 10000 to each codeword in C'). This is a local 
operation, so the code parameters do not change. 



B. Mapping from AC06 to the standard form of CWS 

Three distinct steps may be identified, in building a mapping 
between the AC06 (Af, /) code, and the CWS (Q,C) code in 
standard form, 



(Af,f)^(S A ,C) 



LC 



(G. 



,C')^{Q,C). 



(17) 



First, (Af, /) is re-written as a stabilizer S A and a classical 
code C, using standard definitions. The subscript A on S A 
reminds us that the stabilizer is generated by the generators 
9 A = (91, ■ ■ ■ j 9n), where each generator <j% corresponds to a 
row of Af. Second, a (non-unique) local Clifford transform 
L turns Sa into Qa, leaving C' invariant. Qa is a graph state 
with generators LqaL\ Third, careful choice of appropriate 
generators turn the classical code C into the C used in the 
CWS construction. A fourth issue that arises is the limitation 
on / needed to allow degenerate codes to be considered. These 
three steps and the degeneracy issue are discussed below, one 
at a time. 

1) (Af,f) — ^> (Sa,C): We have already accomplished 
this step by way of Theorem 4, but we review it quickly to 
show the entire chain of steps to achieve standard form. The 
n x In matrix Af describes the generators of a quantum sta- 
bilizer state, which we may denote as Sa, when the left n x n 
half is interpreted as describing X Pauli terms, and the right 
half, Z Pauli terms, following the standard prescription[17]. 
Let the generators of this stabilizer be gA = (91, ■ ■ ■ , g n ); 
each generator gf. corresponds to a row of Af. Let \S) be the 
quantum state stabilized by Sa- 

The Boolean function / defines a classical code, through its 
action on the K bit strings c'j = ji . . . j n ; explicitly, we may 
define 



C' = {c'.|/(c') = l} J 



(18) 



where c'j denotes the complement of (needed because of 
how / is defined in AC06, see Example 2). 

In the CWS standard form, the all-zeros codeword is in the 
classical code C, i.e. the state \S) is in the code. This can be 
arranged by choosing one of the states 1 5 C '. ) in the code and 
applying to the whole code the local Pauli operation that maps 
|5 C '.) to \S). Since this has no effect on the stabilizer Sa, and 
the resulting code is locally equivalent to the original code, 
we now assume without loss of generality that C contains the 
all-zeros codeword. 

2) (Sa,C) — ► (Qa,C): The second step needed is an 
intermediate, but simple map, transforming Sa into graph 
state form[14]. This can be done using Clifford operations on 
individual qubits ("LC transformations"). Importantly, though, 
we must also keep track of how C transforms when the 
stabilizer Sa is transformed, since C is partially defined in 
terms of Sa- 

n 

Let L = Li be the n-qubit operation given by the 

i=l 

tensor product of single qubit Clifford operations Li. When 
transformed by L, the generators of the stabilizer Sa map to 
become 

(gi,...,g n ) -> (g'i,--,g'n) , (19) 

where g'. = LgJJ . Since L also transforms Wj to w'j = 
LwjL\ it follows that the commutation relations of w'^ with g' k 
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are the same as between Wj and g^. Thus, LC transformations 
leave C unchanged, mapping (Sa,C) into (Qa, C). Again, 
just as for Sa, the subscript A on Qa reminds us that the 
generator of this graph state is LgAL\ and originates from 

3) (Qa,C) — ► (Q,C): The final step in transforming the 
quantum code into CWS form involves nailing down a degree 
of freedom which allows C to be changed, without changing 
the stabilizer, or the quantum code specified. In particular, C 
is dependent on the choice of generators for Qa- Let R be a 
binary valued, invertible n x n matrix Rji, which transforms 
a generator set (.91,32, •• -,9n) into (g[,g' 2 , . . .,g' n ), where 

n 

9i H.'/f • (20) 

i=i 

We may keep track of this transform by rewriting Qa as Q, 
though, of course, the stabilizer (and thus the corresponding 
graph) must be left unchanged when the generator set is 
changed. Upon this transformation by R, the code C must 
also be transformed, to keep the quantum code invariant. 
Specifically, if C is written as a K x n matrix, then: 

Theorem 5: The quantum code (Qa, C) is the same as the 
quantum code (Q, C'R). That is, if the stabilizer generators are 
changed by R, the code must also be transformed by matrix 
multiplication by R. 

Proof: We have WjgkWj = (— \Y h gk, and we want to 
calculate j' k given by Wjg' k Wj = (—iy k g' k . Note 

n n 

Wrf'^ = Wj \\ g R kt Wj = J! w 3 g Rkt w 3 

k=l k=l 

n n 

= n(w j9k wj) R *< = n((-iy»g k ) Rkt 

k=l k=l 
n n n 

k=l k=l k=l 

n 

= ((-1)^=^)11^ = (- i y' k 9 , k , 
fc=i 

which gives j' k = ® k=l jkRkt- ■ 

Essentially, this equivalence indicates that row reductions 
in the symplectic n x 2n form of the stabilizer can leave the 
quantum code invariant, if the same row reduction is done to 
the binary code. Moreover, LC equivalence and the choice of 
generators of the graph state do not change the error correcting 
property of the quantum code. Thus, using a row reduction 
transform R, and letting C = C'R, we conclude that (Q,C) 
is a CWS code with dimension and distance identical to the 
original AC06 code (Af,f). 

It must be noted that the row reduction does change the 
errors (in terms of binary strings) detected by the classical 
code. More precisely, for a CWS code (Q,C) in the standard 
form that we have obtained from an AC06 code (Af, /), we 
may define a corresponding (A'j, , /') in the language of AC06, 
by 

f'(cj) = 1, Vcj-eC (21) 
A' = [I A], (22) 



where I is the n x n identity matrix, and A is the adjacency 
matrix of the graph Q. 

The complementary set Csetf of the Boolean function /' 
is no longer the same as the the complementary set Csetf 
of the Boolean function /, but they have same size due to 
the linearity of the transform relating C and C. Moreover, 
given quantum code distance d, the set of induced classical 
error strings Clg(£) for (Q,C) is indeed the AC06 error set, 
specified as {xi,x% ■ ■ -X2k} * w T in Theorem 2 of [12], a 
subset of the complementary set Csetf of /'. 

4) Degenerate codes: The AC06 framework does not dis- 
cuss how to allow for degenerate quantum codes, whereas 
the CWS construction includes these explicitly. The above 
mapping of AC06 to the standard form CWS codes applies 
only to non-degenerate codes, but the method indicates how 
degenerate codes can also be constructed using the AC06 
framework, as follows. Specifically, one must appropriately 
constrain the Boolean function / (ie C). 

All degenerate quantum codes can be expressed using a 
certain form for C , illustrated by the following. Consider a 
degenerate code of distance d, given stabilizer S. Define the 
set 

S d = {E\E e S and wt(E) < d} 

U {-E\E 6 -S and wt(£) < d} , (23) 

where wt(E) gives the weight of the Pauli operator E. If the 
rank of Sd is r, then r independent elements g±,. ..g r £ Sd 
can be chosen, such that (g\, ... , g r , g r +i, ■ ■ ■ , g n ) generate S, 
but g r +i, ■ ■ ■ g n are n °t in Sd- According to the CWS construc- 
tion described in the first step above, these generators imply 
a representation of a classical code C with each codeword 
being for the first r coordinates. In other words, (gi, . . . , g r ) 
stabilizes (Af,f). Due to the one-to-one correspondence 
between / and C, this gives a structure for the values of /, 
from which a search for degenerate codes can initiate. 

C. The algorithm & complexity 

Given the equivalence between AC06 and CWS codes, it 
is insightful to compare the algorithms implied by each for 
finding new codes. Both approaches construct a quantum code 
(Q,C), but each analyze and calculate from different starting 
points. The search algorithm based on the CWS construction 
starts from the analysis of the structure of a given Q, takes 
a specification the desired properties of C, and searches for 
a satisfactory C, eg using the maximum clique algorithm. In 
contrast, the search algorithm based on the AC06 framework 
starts from the analysis of the structure of a given / (ie, C), 
and searches for a stabilizer state Af which is LC equivalent 
to some graph state Q. This is why the two methods are in a 
sense, the mirror image of each other. 

How do the computational complexities of the two ap- 
proaches compare? AC06 implies an algorithm starting from a 
given classical code / to find the quantum code (Af, f). This 
suggests a need to consider 2 2 different Boolean functions. 
In contrast, the cws-maxclique algorithm starts from 2W 
possible graphs (or ideally, a smaller set of just the different 
ones). 
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However, this comparison is incomplete. In practice, if we 
really want to find an particular ((n,K,d)) code, then there 
will be ( 2 K ) classical codes to look at, and for each code 
the AC06 algorithm needs to search for ~ 2 2n possible sets 
of strings. For a given classical code, to check whether a 
particular string is in the complementary set Csetf of the 
code takes K 2 steps. And to check whether a chosen set 
of 2n strings gives a valid stabilizer state [AB] needs n 2 
steps. Therefore, with the AC06 algorithm, the complexity of 
searching for an ((n,K,d)) code is roughly 

n 2 K 2 2 2n2 Qjj . (24) 

This is comparable but slightly worse than the result obtained 
for the CWS-MAXCLIQUE algorithm, in Eq. (7). 

Some simplifications used in CWS-MAXCLIQUE may also 
apply to AC06; in particular, a reduction of the code search 
space due to LC invariance should be considered. In practice, 
in order to find all quantum codes (Af,f), we only need 
to consider the codes C equivalent under column reductions. 
For K > n, this LC equivalence is the same as equivalence 
classification of all the ((K, n')) binary linear codes, where 
n' < n. For fixed n', the number of such codes is given 
by the Gaussian binomial factor ( ,)„ [181. Note 

J V n' ) Gaussian L J 

this classification gives not only all the ((n',K)) codes C 
we need to start with, but also all the ((n',K' < K)) 
codes C . For instance, the ((K = 4, n' = 3)) code 
{(0, 0, 0, 0), (0, 0, 0, 1), (0, 0, 1, 0)}, viewed by column, is an 
((n' =3, if' = 3)) code {(0, 0, 0), (0, 0, 1), (0, 1, 0)}, but not 
an {{n 1 = 3, K = 4)) code. 

IV. The structure theorems 

The ability to search for CWS codes through solving the 
MAXCLIQUE problem is unsurprising; any unstructured search 
problem can be reduced to an NP-complete problem. Thus, 
as it stands, the CWS-MAXCLIQUE algorithm presented in 
Section II is unsatisfactory (at least, for large cases), for the 
search space grows exponentially with the problem size n. 
Moreover, as shown in Section III, the complexity of the AC06 
algorithm is comparably bad, and is thus also unsatisfactory. 

Since a major goal of the study of nonadditive codes is 
identification of codes with parameters superior to all possible 
additive codes, pruning the search space is worthwhile as a first 
step, before applying such brute-force search. 

Is there hope? All nonadditive quantum codes with good 
parameters constructed so far have been CWS codes, as was 
shown in [1]. Also, very recently the ((10,24,3)) CWS code 
was enumerated[3]; this code saturates the linear programing 
bound on code parameters. It thus seems that we should be 
optimistic about finding more CWS codes that outperform 
additive codes. We call an ((n,K,d)) additive quantum code 
optimal if there does not exist any ((n, 2K, d)) additive 
quantum code. One might hope that improved codes could be 
built from optimal ((n,K,d)) additive codes, using the idea 
that these codes could be subcodes of larger (non-additive) 
CWS codes with superior parameters. If this were true, then a 
promising strategy would be to start with the optimal additive 
codes and try to increase the dimension. 



This strategy leads to useful knowledge about the structural 
properties of CWS codes and reveals relations between codes 
with parameters ((n,K,d)) and ((n, K 1 ,d)), where K' > K. 
These relations are especially interesting when given extra 
knowledge about the nature of the classical code C em- 
ployed in the construction. Surprisingly, we find that the low- 
dimensional CWS codes are actually additive. In particular, 
we find that all ((n, 3, d)) CWS codes are subcodes of some 
((n, 4, d)) additive codes. Furthermore, we find restrictions on 
how optimal additive codes can and cannot be subcodes of 
larger CWS codes. 

Before presenting these structure theorems, we review the 
relationship between the linearity of C and the additivity of 

Q = (G,C). 

A. Linearity of C and additivity of Q = (G,C) 

Recall from Theorems 4 and 5 in [ 1 ] that the following facts 
are true: 

Fact 1: If C is a linear code (or equivalently, the word 
operators form a group), then Q = (G,C) is an additive code. 

Fact 2: If Q is an additive code, then there exists a linear 
code C and a graph Q, such that Q = {G,C). 

However, when C is nonlinear, the question of whether 
(G,C) is additive or not is completely open, since it may or 
may not be possible that (Q, C) is local unitary (LU) equivalent 
to some additive code. 

The following example explicitly illustrates this possibility, 
by presenting two CWS codes: (G,C2) with nonlinear C2, and 
(G,Ci) with linear C\. The two codes are LU equivalent to 
each other: 

Example 3: Let 

G = (XZZZ,ZXII,ZIXI,ZIIX) (25) 
Ci = {0000,0110,0101,0011} (26) 
C 2 = {0000,0110,0101,1011}. (27) 

Note that (G,C\) is an additive code since the codewords of 
C\ form a group under binary addition (it is thus a linear 
code). In contrast, since C2 is nonlinear (its set of codewords 
are not closed under addition), (G,C2) is not LC equivalent 
to any additive code. Nevertheless, we can show that Q\ = 
(G, C2) is LU equivalent to Q2 = (G, C\), by giving an explicit 
LU equivalence between the projectors into the two quantum 
code spaces, P\ and ¥2- For this purpose, it is convenient 
to first transform by H234 = H2 (E> H3 ® H4 and disregard 
normalization factors, such that 

P[ = H234P1H234 

= I + XXXX + YYYY + ZZZZ (28) 

P-2 = H234P2H234 

= I + ZZZZ 

+ ^(XXXX + YYYY + XXYY + YYXX 
-XYYX - YXXY - XYXY - YXYX) (29) 

From Theorem 4.2 of [19], LU equivalence need only consider 

U = Ui <g> U 2 ® U3 ® U4 where Ui maps X to aX + bY and 
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Y to bX - aY. We find that UP[W = P£, if U is denned 
such that 

C/^[/t = [Xi- (-l)^Yi\/V2 (30) 
C/.y.C// = [^ + (-l)LV2Jy.]/V2, (31) 

where [*/2j is for i < 2 and 1 otherwise. The existence of 
this LU equivalence is unsurprising, since it is known [20] that 
any ((4, 4, 2)) code is LU equivalent to the additive [[4, 2, 2]] 
code. 

In general, for a CWS code Q = (Q,C) with a nonlinear C, 
we cannot directly infer that Q is nonadditive. However, for 
fixed n and d, if we seek a code with optimal K and only find 
((rt, K' > K, d)) codes Q = (G,C) with nonlinear C, then we 
can conclude that Q nonadditive. Put another way, if we fix n 
and d, do an exhaustive search over all the graphs and classical 
codes, and only find quantum codes with nonlinear classical 
codes C for the optimal ((n, K, d)) CWS codes, then we can 
conclude that the optimal ((n,K,d)) CWS codes we found 
are indeed nonadditive. This can be shown by contradiction: 
if Q = (Q,C) is additive, then there exists some local unitary 
operation U — ®"=i where each Ui is a single qubit 
operation, such that UQU^ = Q' and Q' is additive. Then, 
according to Fact 2, there exists a linear code C and a graph 
Q' such that Q' = (Q',C). 

B. Structure theorems 

We now present and prove some structure theorems govern- 
ing CWS codes, and provide several useful corollaries. Recall 
that we say an additive ((n, K, d)) quantum code is optimal 
if there is no ((n, 2K, d)) additive quantum code. 

Our first theorem concerns CWS codes with dimension 2: 

Theorem 6: All ((n, 2, d)) CWS codes are additive. 
Proof: By the CWS construction, an ((n,2,d)) CWS 
code is spanned by basis vectors of the form {wi\S), W2\S)}, 
with word operators w\ = I = Z Cl , w 2 = Z C2 . However 
{wi, W2} form a group. So according to Theorem 5 of [1] (or 
Fact 1), this CWS code is an additive code. ■ 

A natural corollary of Theorem 6 is 

Corollary 1: If an additive code of parameters ((n, 1, d)) 
is optimal, then there do not exist any CWS codes with 
parameters ((n, K > l,d)). 

From corollary 1, it follows that the ((7, 2, 3)) and ((9, 2, 3)) 
nonadditive codes given in [21] and the ((11, 2, 3)) code given 
in [19] are not local unitary (LU) equivalent to any CWS 
code, for they are not LU equivalent to any additive code. 
This implies that there exist codes that are outside the CWS 
construction, as was claimed in Fig. 1. 

Now we present a theorem concerning CWS codes of 
dimension 3: 

Theorem 7: Any ((n, 3, d)) CWS code is a subcode of some 
((n, 4,d)) stabilizer code. 

Proof: By the CWS construction, any ((n,3,d)) CWS 
code has the form (Q, C\) with C\ = {ci =0, c 2 , C3}. Consider 
a new code (<?,C 2 ) with C 2 = {ci =0, c 2 , c 3 , c 2 © c 3 }. From 
Theorem 1, it follows that C\ detects errors in Clg(£). To 
prove Theorem 7, we need to show that C 2 also detects those 
errors. It is clear that C 2 is a group with generators c 2 , C3 and 



that c 2 ©c.3 ^ Clg{£) because c 2 ffi(c 2 ©C3) = C3. Therefore 
C 2 detects all of Clg{8). Theorem 1 also requires that for each 
E e £ either Clg(E) ^ or for all i, Z Ci commutes with 
E. The latter constraint is satisfied by C 2 since Z C2 ® C3 E = 
Z C2 Z Ca E = EZ C2 Z C3 . Finally, since {I, Z C2 , Z c \ Z C2ffiC3 } 
is a group (and thus a linear code), according to Theorem 5 
in [1] (or Fact 1), this CWS code is a stabilizer code. ■ 

Two natural corollaries of Theorem 7 are: 

Corollary 2: If an additive code of parameters ((n, 2,d)) 
is optimal, then there do not exist any CWS codes with 
parameters ((n, K>2, d)). 

Corollary 3: There does not exist any ((7, 3, 3)) CWS code, 
even though the linear programing bound does not rule out this 
possibility. 

The two structure theorems above imply that CWS codes 
with parameters better than the optimal ((n,K,d)) additive 
codes need dimension K > 4. We do know examples where 
K = 4, as the ((5, 6, 2)) code [5] and the ((5, 5, 2)) code [6] 
beat the optimal additive code with parameters ((5, 4, 2)) [22]. 

Theorem 7 says that a CWS code of dimension 3 is a 
subcode of some additive code with higher dimension. This 
invites a related question: when might an optimal additive 
code, of dimension K, be a subcode of some CWS code of 
higher dimension? Unfortunately, we can show that in some 
sense, optimal additive codes cannot be subcodes of larger 
CWS codes, though we cannot show the impossibility in the 
most general setting, due to the fact that C may be nonlinear 
even if a CWS code is additive. 

Motivated by LU equivalences like the one demonstrated in 
Example 3, we show that if C\ is a linear code, then an optimal 
additive code (Q,Ci) cannot be a subcode of any CWS code 
(G,C 2 ), where Ci C C 2 : 

Theorem 8: Given a CWS code (<7,Ci) with parameters 
((n, K, d)), if B is a linear subcode of C containing J < K 
codewords, then there exists an additive code (<5,C 2 ) with 
parameters ((n, K' = 2J,d)). 

Proof: By the CWS construction the classical codewords 
C\ = {ci, c 2 , . . . ca'} of (G,Ci) can be arranged such that 
ci = 0. From B construct the linear classical code C 2 = 
{bi,b 2 ...bj, v ffi bi, v ffi b 2 .. .v © bj} where v e C x but 
v ^ B. Then (5,C 2 ) is clearly an n-qubit CWS code with 2 J 
codewords. It is an additive (stabilizer) code by Theorem 5 of 
[1] since C 2 is a group. 

It remains to check the error-correction conditions. Theo- 
rem 1 ensures that C\ detects errors in Clg(£), i.e. no error 
can turn one codeword into another: 

Ci © Cj © e ^ for all e £ Clg{£) . (32) 

The same condition for C 2 is 

b, © v k © bj © v' © e ^ , (33) 

where k, I £ {0, 1}. Since the 6s are a group this reduces to 

bi © v fe © e ^ (34) 

which is true, due to Eq.(32), and the fact that bj,0,v € C\ 
for all i. 

Theorem 1 also tells us that for all E € £ either (a) 

Clg {£) ^ or (b) for all i, [Z c - , E] = 0. (Q, C 2 ) has the same 
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graph Q as (G,Ci) so whenever (a) is satisfied for (G,Ci) it 
will be for {0,0%). For C2 (b) becomes for all i = 1, J and 
£• = 0,1 [Z b 'Z vfc ,£:] = 0. Again, since b;,v e C x for all i, 
this is condition is met. ■ 

Corollary 4: An optimal additive code (G,C) (for which C 
must be linear) cannot be extended to become a larger CWS 
code merely by adding codewords to C. 

Proof: If the code could be extended in this way, by 
adding even just one vector, then there would exist an additive 
code with twice as many vectors and the same distance as the 
original code. This contradicts the statement that the original 
code is optimal. ■ 

These structure theorems rule out certain strategies for 
finding non-additive codes with parameters superior to addi- 
tive codes, but suggest other approaches. Since an additive 
((n,K,d)) code (G,Ci) must have linear C\, Theorem 8 and 
corollary 4 tell us that in practice we cannot search for an 
((n, K' > K, d)) CWS code (G,C 2 ) just by adding codewords 
to C\. However, Example 3 hints that we may be able to 
shoehorn an optimal ((n,K,d)) additive code into a CWS 
code (G,C) with nonlinear C, via some LU transform. This 
gives hope to a strategy of adding codewords to C to search 
for ((n, K' > K, d)) CWS codes; such hope suggests that it is 
worthwhile both to further explore conditions under which two 
CWS codes can be linked by an LU transform, and to better 
understand the structural properties of CWS codes constructed 
from nonlinear codes. 

V. Discussion 

CWS-MAXCLIQUE is an algorithm which may be usefully 
employed in the search for new quantum codes, both additive 
and non-additive, as described by the CWS construction. 
Given n and K, the algorithm can be used to search for 

an ((n,K,d)) code (G,C), with a complexity which grows 

2 

roughly as 2" . In practice, by employing a number of search 
space simplifications, by pruning the set of graphs G to explore 
based on LC equivalences, and by taking guidance from 
structural theorems about CWS codes, CWS-MAXCLIQUE and 
randomized variants of it have been used realistically[l] to 
explore codes with parameters up to n = 11 and K = 32. 

Many interesting questions arise in the construction of this 
algorithm. For example, it is likely that CWS-MAXCLIQUE 
can be improved with more memory efficient implementa- 
tions; reductions to other NP-complete problems may also 
allow faster exploration of specific search spaces. Moreover, 
many of the simplifications used in CWS-MAXCLIQUE should 
also be applicable to the algorithm introduced by the AC06 
framework; and in return, any code isomorphisms useful in 
simplifying AC06 should apply to CWS-MAXCLIQUE. 

CWS codes present a rich structure, only partially described 
by the three structural theorems presented here. We believe 
that there are promising strategies for identifying new non- 
additive quantum codes based on expanding known additive 
codes, but such a strategy has to be executed carefully, 
because of limitations imposed by the theorems. Nevertheless, 
given an optimal ((n,K,d)) additive code, there is hope for 
success with a strategy of adding codewords to C to search 



for ((n,K' > K,d)) CWS codes, because of potential LU 
equivalences with some non-additive code. This hope suggests 
that it is worthwhile both to further explore conditions under 
which two CWS codes can be linked by an LU transform, and 
to better understand the structural properties of CWS codes 
constructed from nonlinear codes, so that more new quantum 
codes can be found. Indeed, one successful application of 
this idea results in new CWS codes encoding several more 
qubits than the best known codes [4]. It is an open question 
to determine if these nonadditive "quantum Goethals-Preparata 
codes" are LU equivalent to any additive quantum code. 

Finally, despite the encompassing success of the CWS 
construction in describing all known non-additive codes with 
good parameters, we point out that there do exist codes, such 
as ((7,2,3)) and ((9,2,3)) codes, which are outside of the 
CWS construction. Since these codes are not LU equivalent 
to any CWS code, further new ideas will need to be developed 
to reach outside the stabilizer framework, for a complete 
understanding of quantum error correction codes. 
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